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Gamma-Ray Bursts (GRBs) are now considered as relativistic jets. We analyze the gravitational 
waves from the acceleration stage of the GRB jets. We show that (i) the point mass approximation 
is not appropriate if the opening half-angle of the jet is larger than the inverse of the Lorentz factor 
of the jet, (ii) the gravitational waveform has many step function like jumps, and (iii) the practical 
DECIGO and BBO may detect such an event if the GRBs occur in the Local Group. We found that 
the light curve of GRBs and the gravitational waveform are anti-correlated so that the detection of 
the gravitational wave is indispensable to determine the structure of GRB jets. 



I. INTRODUCTION 

The gravitational wave astronomy enters its first step. The ground-based laser interferometers such as TAMA300, 
LIGO and GEO600 have already begun operations covering 10 Hz - kHz band, and VIRGO will be also in operation 
soon. The Laser Interferometer Space Antenna (LISA) [lj covering 10 -4 - 10 -2 Hz is expected to be launched in 
2011. Recently, in addition, the Decihertz Interferometer Gravitational wave Observatory (DECIGO), which aims the 
frequency band around decihertz (10 -2 - 10 Hz), has been proposed [2|. The decihertz antenna is also considered in 
the NASA Roadmap Beyond Einstein and is called the Big Bang Observer (BBO) 0). 

It might be possible to construct a practical DECIGO and BBO with h lms ~3x 10~ 24 at ~ 0.1 Hz. Although this 
limiting sensitivity is ~ 10 3 times less than the ultimate quantum limit h rms ~ 3 x 10~ 27 at ~ 0.1 Hz for a 100 kg 
mirror mass, the practical DECIGO and BBO may achieve the following important results (1) The chirp signals 

of coalescing binary neutron stars at z ~ 1 may be detected with a signal-to-noise ratio (SNR) ~ 10 for one year 
observation. (2) The time variation of the Hubble parameters of our universe may be determined by the chirp signals 
so that we can directly measure the acceleration of the expansion of the universe. (3) The stochastic gravitational 
waves from inflation may be detected if flow ^ 10~ 15 , which is the upper limit from cosmic microwave background. 
(4) An arcminute position and an accurate (~ 0.1 s) coalescence time of the coalescing binary neutron stars at ~ 300 
Mpc may be determined a week before the final merging, so that the all band electromagnetic detectors can perform 
the simultaneous follow-up observation. An accurate position may be also important for the reconstruction of the 
density distribution and the shape of the galactic dark halo 0, la • 

In this paper, we consider gamma-ray bursts (GRBs) as possible sources of the gravitational waves for DECIGO 
and BBO. The GRBs are the most luminous objects in the universe. The recent discovery of the afterglow has 
revolutionized the understanding of the long duration GRBs (e.g. 0>I3)- We are now confident that long GRBs with 
~10s duration eject relativistic flows with Lorentz factors 7 > 100 and the ejecta is collimated in a jet. This means 
that the nonspherical acceleration of a lar ge m ass is associated with GRBs. Therefore it is guaranteed that the GRB 
jet acceleration emits gravitational waves [Tol UlI ll2T| . We will show that the typical frequency of the gravitational 
waves is about the inverse of the duration of GRBs ~ 0.1 Hz, so that the decihertz band of DECIGO and BBO best 
suits the detection of gravitational waves from the GRB jet f§ llllf) . 

Here we concentrate on the gravitational waves from the acceleration stage of GRB jets, while we also expect 
the gravitational waves from the formation of the central engine of GRBs [R3, [lj, [ija, [Tg that could be detected 
by the ground-based laser interferometers. We know that the accelerated jet exists although we do not know how 
it is accelerated and what structure it has. Since the jet is accelerated near the central engine where the optical 
depth is quite large, the electromagnetic probe is useless. The high-energy neutrinos from GRBs (\17\) may provide 
a indirect probe of the central engine, while the gravitational waves may bring us the direct information of it. Thus 
the gravitational waves could be an ultimate way to see the acceleration regime. 

The gravitational waves from the GRB jet have a "memory" , that is, the metric perturbation does not return 
to its original value at the end of the jet acceleration [l^. The gravitational wave memory is the change in the 
transverse-traceless (TT) part of the Coulomb-type (oc 1/r) gravitational field generated by the masses of the system. 
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These masses could be stars, neutrinos |l9t l20j . gravitons and in our case the GRB jets. 

Segalis and Ori [lOj analyzed the gravitational wave emitted when an ultrarelativistic blob is ejected from a massive 
object by using a point mass approximation. They showed that the gravitational wave emission from such a source 
is anti-beamed. This contrasts with the electromagnetic emission, which is beamed into the forward direction. They 
also obtained the detectable distance of about 15 Mpc and estimated the detection rate of about 1 per year by the 
advanced LIGO, assuming that the total energy of the GRBs is about ~ 10 52 erg and that GRBs follow supernovae. 

On the other hand, recent studies suggest that the jet structure is essentially important for understanding the 
GRB phenomena |2J, |25j, |26( . Depending on the viewing angle, the structured jet model might be able to explain the 
diverse phenomena of GRBs, such as short GRBs, long GRBs, X-ray flashes (XRFs) and X-ray rich GRBs [27|. Here 
XRFs and X-ray rich GRBs are fast X-ray transients that are considered to be related to long GRBs [28l l29l ] . Little is 
known about short GRBs since no afterglow was observed although many suggestions to show the similarity to long 
GRBs exist. 

It is difficult to determine the jet structure only by the electromagnetic waves, since the photons are beamed into 
the forward direction and we can observe only a small angular portion ~ 7 -1 ~ 0.01 rad of the whole jet. On the 
other hand, since the gravitational wave emission is anti-beamed, it can be observed at larger viewing angles, so that 
we can see wider angular portion of the whole jet. Therefore the information carried by the gravitational waves is 
important to know the jet structure. In this paper we investigate the gravitational wave from an inhomogeneous jet. 
We model the inhomogeneous jet in which multiple sub-jets with a finite opening half-angle are ejected to various 
directions. The jet of GRBs should not be a single homogeneous one since a single explosion cannot produce a variable 
GRB. We discuss what kind of information is carried by the gravitational wave observed by DECIGO and BBO. We 
also discuss the detectability of such events. 

This paper is organized as follows. In Sec[HJ we consider the gravitational wave from a single jet acceleration. We 
investigate the angular dependence of the gravitational wave from a jet with a finite opening half-angle. In Sec lIIII 
we discuss the gravitational waveform from the GRB jet adopting a unified model of the GRB proposed by some of 
authors j^. In Sec. IIVI we calculate the strain amplitude detected by the gravitational wave detector. We also 
evaluate the maximum detectable distance to the GRBs for such events by using DECIGO and BBO. Finally, we 
summarize this paper and discuss our results in Sec. We use the geometric units G = c = 1 and the signature 

(- + ++)■ 



II. GRAVITATIONAL WAVE MEMORY OF A SINGLE JET 



In this section we consider the gravitational waves from a single jet acceleration. Actually a single explosion cannot 
produce a variable GRB. Many jets are to be ejected since the total duration of the burst is usually much longer than 
the variability timescale jsjj. Many jets case will be discussed in the next section. 

First let us remember the generic form of the gravitational wave memory ^[ [2^ • We consider that the system 
consists of freely moving masses before and after the burst. Then the variation of the metric perturbation (which is 
the gravitational wave memory) is given by 

Here A means the difference between before and after the burst, wia is the mass labeled by A(= 1, 2, • • •), r is the 
distance from the source to the observer, (3 A is the velocity of mass A, 7,4 = (1 — Pa)^ 1 ^ 2 is the Lorentz factor of mass 
A, and 9a is the angle between (3\ and the direction from the source to the observer. Small Latin indices ■ ■ ■) 
run over spatial coordinates x, y, z, and the quantities r, f3 A , 9a are measured in the lab frame. The TT part of Ah%j 
is obtained by AhJ? = P ik PjiAhki - (l/2)P ij (P M Ah k i) where Pij — $ij ~ UiUj is the projection operator and n.$ is 



the unit vector from the source to the observer |3JJ| . The source quantities are to be evaluated at the retarded time. 

The gravitational wave memory of a point mass acceleration is easily calculated from equation (|2.1(l lTl UlI IT3 | . Let 
us consider that a point mass M is accelerated to a Lorentz factor 7. Conservation of energy requires that the mass 
of the accelerated particle is m — M/7. We take a coordinate such that the z axis is the direction from the source to 
the observer and (3 % oc (sin 9 cos (j), sin 9 sin 0, cos 9) . Then equation j|2.1|l gives the amplitude of the gravitational wave 
memory as 

27771/3 2 sin 2 9 , 

Ah = Ah++iAh x = ■= ^e 24 , 2.2 

r 1 — p cos 9 

where Ah + = AhJJ = —Ahyy and Ah x = Ahjy = AhyJ are two polarizations of the gravitational waves, plus and 
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cross modes. For 7>1 and 9 <C 1, we have 



\Ah\ 



(2.3) 



while for 7 > 1 and > 7 1 



\Ah\ 



27m 



(1 + cos 0). 



(2.4) 



Therefore the gravitational wave memory is anti-beamed with a similar amplitude in almost all direction but little 
amplitude in the forward direction 9 < 7" 1 - The waveform has a form like a step function. As we can see from 
equation (12. 2[) , the rise time of the amplitude from its initial to final value is about the time to reach (3 ~ 1 if the 
energy 7m is conserved. Thus the rise time is about St ~ 10 7 cm/c ~ 10~ 3 s in the fireball model of the GRBs 
[32I [33 • If a counter-jet exists [34| , we sum both contributions to obtain 



|AA| 



47m 



1.6 x 10 



-22 



27m 



3 x 10 51 erg J V IMpc 



(2.5) 



for 7 ^ 1 and 9 ^> 7 _1 . (In this paper, however, we do not consider the counter jet |52l|.) 

In reality the opening half-angle of the GRB jet is larger than 7 -1 [37l l3S| | . Thus equation l|2.2|l is not appropriate 
for the GRB jet. To examine the effects of the finite opening half-angle, let us first consider a simple model, that is, 
an uniform thin jet with an opening half-angle A9 that is accelerated to a Lorentz factor 7. Since linear perturbation 
is additive, the gravitational wave memory of the GRB jet is given by 



Ah = I s[nd ^H(A9 



)H 



COS( 



cos A9 — cos 9 V cos 9 
sin 9„ sin 9 



Ah., 



(2.6) 



where Ah. is the amplitude of a point mass in equation (|2.2|l . H (x) is the Heaviside step function describing that 
the jet is inside the cone, 9 V is the viewing angle that the jet axis (in the xz plane) makes with the z axis, and 
' "' ' " ' " "' ' 2n(l — cos A9) ~ 7r(A6*) 2 is the solid angle of the jet. Replacing <p with — in equation 



(2.7) 



An = Jo 71 Jo sin 

(12. 6|) . we can show A/i_ = 0. For Ah + we can simplify equation 1)2. 6fl as 

~ Ae 27m/? 2 sin 3 9 sin 2A<^> 



where 



Ah+ = 



A<j> = 



d9- 



max[O,0„— A6] 



rAtt 1 — (3 COS ( 



: AO > 9 V and < 



iA8-c 



<A9-9, 
: others. 



(2. 



As in the point mass case, the gravitational waveform becomes like a step function, and the rise time of the amplitude 
is about the time to reach (3 ~ 1 if the energy 7m is conserved. 

In Figure ^ the amplitude of the gravitational wave memory for the GRB jet in equation (|2.7|) is shown as a 
function of the viewing angle 9 V by solid lines. We adopt the Lorentz factor 7 = 100. The opening half-angle of the 
jet is AO — 0.01,0.1 and 0.2 from top to bottom. We can see that, because of the anti-beaming effect, the amplitude 
is suppressed in the forward direction, that is, when the viewing angle is smaller than the opening half-angle of the 
jet 9 V < A9. The amplitude of the memory for the point mass in equation (|2.2(l is also shown by the dashed line. 
We can find that the point mass is a reasonable approximation of the jet when the opening half-angle is smaller than 
the inverse of the Lorentz factor A9 < r )~ 1 or the viewin g an gle is larger than the opening half-angle 9 V > A9. The 
opening half-angle of the jet is usually small A9 ~ 0.1 [37J, |38( so that we observe the jet from off-axis direction 
usually. This means that the effects of the finite opening angle may not reduce the event rate of the gravitational 
waves from GRBs so much. 



III. GRAVITATIONAL WAVEFORM FROM GAMMA-RAY BURST JETS 



A. Typical frequency 



The leading model of the GRB emission is the internal shock model 30] . The internal shock occurs in the relativistic 
wind when the fast moving flow catches up the slow one. The wind can be modeled by a succession of relativistic 
shells. A collision of two shells produces a single pulse light curve, whose superposition makes a whole light curve. 



4 




FIG. 1: The amplitude of the gravitational wave memory for the GRBs jet in equation l|2.7|l is shown as a function of the 
viewing angle 8 V by solid lines. We adopt the Lorentz factor 7 = 100. The opening half-angle of the jet is A8 = 0.01, 0.1 and 
0.2 from top to bottom. We can see that the amplitude is suppressed in the forward direction 8 V < A8. The amplitude for the 
point mass in equation 112.21 is also shown by the dashed line. Clearly the point mass approximation is good when the opening 
half-angle is smaller than the inverse of the Lorentz factor A8 < 7" 1 or the viewing angle is larger than the opening half-angle 
9 V > A8. The amplitude is normalized by h = 7m/? 2 /r ~ 2.7 x 10~ 23 (7m/10 51 erg)(r/lMpc)- 1 . 



When the GRB ejects many jets, the gravitational waveform from the GRB differs from a single step function, but 
has many steps since linear perturbation is additive. There are three timescales that determine the temporal structure 
of the gravitational wave memory of the GRB. The first is the rise time of the gravitational wave amplitude due to a 
single jet AT r i so . As argued in the previous section, AT r ; se is related to the acceleration time of the jet. The second is 
the separation between successive jets AT sep ~ 0.1s. The third is the total duration of the burst AT, which is about 
AT ~ 10 s for long GRBs. Then the typical order between three timescales is 

AT risc « AT sop « AT. (3.1) 

Therefore the typical gravitational waveform should have many steps. Each step corresponds to the ejection of a jet, 
and the final amplitude of the gravitational wave memory is mainly determined by the total energy of the jets ejected 
from the GRB. The final amplitude is reached after the total duration AT. Thus the typical frequency in Fourier 
space of the gravitational waves is expected as 

^ c ~ a~t ~ °- 1 Hz ' (3 - 2) 

This band is just the decihertz band of DECIGO and BBO. 



B. A unified model of GRBs 



Recently some of authors proposed a unified model of the GRBs [2j| . In this model, the central engine of the long 
and short GRBs, XRFs and X-ray rich GRBs is the same, and the apparent differences come essentially from different 
viewing angles. We assume that the jet of the GRB consists of multiple sub-jets. These sub-jets have angular size of 
~ 7 _1 <C 0.1 and are distributed within the whole GRB jet (~ 0.1). If many sub-jets point to our line of sight, the 
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FIG. 2: The angular distribution of AT to t = 350 sub-jets confined in the whole GRB jet in our simulation. The whole jet has 
the opening half-angle of A# tot = 0.2 rad, while the sub-jets have A# su b = 0.02 rad. The axes and the angular size of sub-jets 
are represented by crosses and the dotted circles, respectively. "A" represents the center of the whole jet and is hidden by the 
lines of sub-jets. (Also refer to Yamazaki, Ioka and Nakamura, Astrophys.J. 607 (2004) L103-L106 |23|.) 



event looks like a long GRB, while if a single sub-jet points to us, it looks like a short GRB. If our line of sight is 
off-axis to any sub-jets, the event looks like an XRF or X-ray rich GRB. 

An example of the sub-jet angular distribution is shown in Figure where iVtat = 350 sub-jets are ejected for 
AT = 30 sec and each sub-jet has 7 = 100 and A0 su b = 0.02 rad. When we observe the source from the 8 V = axis 
(case "A"), we see spiky temporal structures (Fig. 3-A) and may identify the case "A" as the long GRBs. When the 
line of sight is away from any sub-jets (cases "Bi" and "B2"), soft and dim prompt emission, i.e. XRFs or X-ray rich 
GRBs are observed (Fig. 3-Bi, 3-B 2 ). If the line of sight is inside an isolated sub-jet (case "C"), its observed pulse 
duration is much smaller than the case "A" (Fig. 3-C), so that we may observe the short GRB. The details about the 
calculations of the X-ray and gamma-ray light curves in Figure are given in 

In Figure^ we show the plus mode of the gravitational wave as a function of time assuming the same model as 
above. Each panel ("A", "Bi", "B2" and "C") corresponds to the viewing angle presented in Figure|5]and the X-ray 
and gamma-ray light curves in Figure [3J From Fig.l for Ad ~ 7 j the point mass approximation is good. Since 
A# su b = 27 _1 in our case we use the point mass approximation for each sub-jet and neglect the rise time of the 
amplitude, for simplicity. Each step corresponds to the ejection of a sub-jet, and the whole gravitational waveform 
has JVtot = 350 steps. Since the sub-jets are ejected for AT = 30 sec, the final amplitude is reached after the duration 
AT. 

We find that the final amplitude in case "A" is much smaller than the other cases and goes up and down more 
frequently than the others. This is because each memory of the sub-jet is proportional to e 2l< ^ from equation 12.21) 
and the distribution of cf> is nearly isotropic for an on-axis observer like case "A" so that each memory tends to be 
canceled. On the other hand, if the viewing angle is large, the distribution of cf> is not isotropic and either of two 
polarizations dominates the other one. Therefore, the amplitudes in cases "Bi","B2" and "C" increase monotonously. 
In addition, their waveforms resemble each other because the gravitational wave emitted by each sub-jet is nearly 
isotropic oc (1 + cos 8) for 9 < 0.2. 

Comparing the gravitational waveform with the gamma-ray one in Figure |3 we can see that the gravitational 
wave amplitude is roughly anti-correlated with the photon flux. In case "A", the gamma-ray is strong while the 
gravitational wave is weak. In cases "Bi" and "B2" , the gamma-ray is weak while the gravitational wave is strong so 
that we may confirm the existence of many off-axis sub-jets which are not identified by the gamma-ray. In case "'C', 
the gamma-ray is short and strong while the gravitational wave is long and strong. This means that the gravitational 



6 




20 30 40 50 25 30 35 

T obs (sec) T obs (sec) 



FIG. 3: The observed X-ray and gamma-ray photon number flux from the multiple sub-jets, corresponding the cases "A" (the 
upper left), "Bi"(the upper right), "B2"(the lower left) and "C"(the lower right) in Figure|5| Here the flux in 30 — 400 keV 
band for each case is shown. (Also refer to Yamazaki, Ioka and Nakamura, Astrophys.J. 607 (2004) L103-L106 27].) 



waves carry information on the off-axis jets. Thus, the gravitational waves are indispensable to explore the structure 
of the jet. For instance, in other models of the XRFs, emission comes from on-axis jets, so that the gravitational wave 
would have less amplitude. In other words, the gravitational wave can be used to verify or refute models of GRBs 
including our unified model of GRBs. 



IV. DETECTABLE DISTANCE 



For the detector with a peak sensitivity h pea k at frequency f pea k , the optimal sensitivity to the gravitational wave 
memory is ~ h pea k as long as the rise time of the memory t m is smaller than ~ l// pe afc [lB> H3- For example, we 
consider a simple waveform as follows: 

( : t < 0, 

h(t) = I Ah m {t/t m ) :0<t< t m , (4.1) 

I A/i m • t > tmi 

which we call the slope waveform. We can calculate the Fourier component of this waveform, h = htye^^dt, 
which yields [&| 

\M 2 = ^^(l-™s2irft m ). (4.2) 

Note that \h\ 2 ~ (A/i m /27r/) 2 for / <C l/t m . In Figure03 we plot the characteristic amplitude of the gravitational 
wave memory, h c (f) — 2f\h\ for t m = 30 sec by the dashed line. We find that h c (f) is constant at / < l/t m , while it 
decreases for / > l/t m . Therefore the detector has the maximal sensitivity /i pea k for the gravitational wave memory 
as long as l/t m < / bs- We also plot h c (f) for the multiple step waveform (case "B 2 ") by the solid line. Here we 
fix the total energy and the distance of the source as 2mj = 3 x 10 51 erg and IMpc, respectively. For this case, the 
amplitude is larger than that for the slope waveform at / > 1/ 1ST. This is because the waveform in the sub-jet model 
has the smaller timescale, AT sep . 

In Figure© we plot the noise amplitude h n (f) = [5fS h (f)}^ 2 for LIGO II, LISA and DECIGO/BBOHijWhere 
Sh(f) is the spectral density of the detector's strain noise at frequency /, and 5 comes from sky-averaging |4l|. We 
also plot the characteristic amplitude of the gravitational wave memory by the bold line. For the low-frequency 
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FIG. 4: The waveform of the gravitational wave memory from the GRB is shown. There are three timescales that determine the 
temporal structure of the waveform, the total duration AT ~ 30 s, the separation between successive jets AT scp ~ AT/N to t = 
0.086 s and the rise time of the amplitude due to a single jet AT acc ~ 10 -3 s (although we neglect the rise time here). The 
final amplitude is reached after the total duration AT. The normalization, ho is the same one in Figure Q 

region, / < 1/AT, we use the result for the slope waveform with Ah m = 4m^/3 2 /r, while we plot it by the bold 
dashed line in the high-frequency region, / > 1/AT, because it depends on the structure of the GRB jet. The upper 
bold line corresponds to a GRB with (AT, r, 2jm) — (30sec, lOkpc, 3 x 10 51 erg), and the lower one to a GRB with 
(AT,r,2jm) = (lOsec, IMpc, 3 x 10 51 erg). The signal-to-noise ratio (SNR) is given by 

sNR2 =r*«w' <4 - 3) 

thereby one can roughly read off the SNR from Figure When the duration of GRBs is shorter than lOsec, the 
detectable distance reaches the Local Group (~ IMpc) by using the practical DECIGO and BBO. On the other hand, 
the detectable distance is reduced when the duration is longer than lOsec. If a GRB happens in our galaxy, we can 
get the information on the high-frequency region / > 1/AT. This may make it possible for us not only to detect the 
gravitational wave memory, but also to extract important information about the structure of the GRB jet from the 
observed gravitational wave. 

V. SUMMARY AND DISCUSSION 

In this paper, we have investigated the gravitational wave emitted in the acceleration phase of the GRB jet. At 
first, we have considered the angular distribution of a GRB jet with a finite opening half-angle Ad. As a result, we 
have shown that we can deal with the jet as a point mass if the opening half-angle is smaller than the inverse of 
the Lorentz factor, 7 _1 , or the viewing angle 9 V is larger than the opening half-angle. We have also shown that the 
gravitational radiation is weak in the angular size ~ 7~ 1 . Next, we have considered the waveform of the gravitational 
wave from the GRBs. We have found that many steps usually appear in the waveform. This has been explicitly 
shown for the multiple sub-jet model of GRBs. If we observe the GRBs from the front, the observed amplitude of the 
gravitational wave is suppressed because the wave from each sub-jet tends to cancel one another. On the other hand, 
the amplitude increases monotonously for the off-axis observer. In addition, we have considered the detectability of 
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FIG. 5: The characteristic amplitude h c of the gravitational wave memory is plotted as a function of frequency for the whole 
GRB jets. The solid line represents h c for the multiple step waveform (case "B2") with the total energy 27m = 3 x 10 51 erg, 
the duration AT = 30 sec and the distance to the source r — IMpc. The dashed line represents h c for the slope waveform 
with the duration t m — 30 sec. Here the amplitude Ah m in equation 14.11 is chosen so that both amplitudes coincide in the 
low- frequency limit. 



GRBs by using the practical DECIGO and BBO. Wc have found that we may detect the gravitational wave emitted 
from the GRB jet if it locates within 1 Mpc. 

Since the gravitational wave memory is anti-beamed, almost all GRBs detected by the gravitational waves will be 
observed with off-axis viewing angle. The photon emission from such off-axis GRBs is softer and dimmer than the 
on-axis ones due to the relativistic beaming 0> 03 ■ Therefore the GRBs detected by the gravitational wave memory 
may be observed as XRFs or rather softer events in the UV/optical bands [34| . The prompt electromagnetic 

emission is too dim to be observed in the case of large viewing angle. Even then, bright radio afterglows can be seen 
several months after the GRB explosion since the relativistic beaming effect becomes weak [3^, |3(| • The association 
with a bright radio source might help us to identify the signal of the gravitational wave memory. In addition, the 
association with supernovae observed in the optical band may be also helpful. XRFs may be GRBs with large viewing 
angles [53, EH ■ The gravitational waves will be observed prior to the XRFs with the time delay being longer for softer 
and dimmer XRFs [42}. Thus the gravitational waves allow one to test the off-axis GRB model for the XRFs. 

Unfortunately the event rate of the gravitational waves associated with the detectable GRBs may be small. Since 
the true GRB rate is estimated as ~ 250 Gpc -3 yr -1 j3?|], the event rate within the detectable distance ~ 1 Mpc (see 
§ |Wj| is only 3 x 10 -7 yr , that is very small compared with the detection rate of the gravitational wave by the central 
engine However the gravitational wave event rate without detectable GRBs could be much higher. Recently, 

many observations support a massive stellar origin for the long GRB 0, [i| and it is widely assumed that a massive 
collapsing star ( "collapsar" ) is the progenitor [47ll48|. Suppose the GRB jets have the same total energy. Although 
the GRB jet cannot break through the star because of weak collimation or it is baryon-loaded, so that 2 < 7 < 30, the 
resultant gravitational wave may have same amplitude as usual succeeded GRBs. Hence the detection rate of such 
events may be improved. Furthermore we calculate the detection rate for mean events. However the total energy is 
distributed around these mean values. Then the detections will be dominated by rare, energetic events. For instance, 
even if only ~ 2% events are 10 times more energetic, the detection rate will be increased by ~ 2% x 10 3 ~ 20, since 
we can observe 10 times farther events. Even if we take account of the above, the detection rate is 6 x 10 -3 yr" 1 
optimistically. We need ~ 10 times better sensitivity than that of practical DECIGO/BBO for the detection rate ~ 1 
yr" 1 . 
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FIG. 6: The noise amplitude h n of LIGO II, LISA and DECIGO/BBO is also plotted by the dotted line. The characteristic 
amplitude for the gravitational wave memory of GRBs is also plotted by the bold solid line. The upper bold solid line 
corresponds to a GRB with the total energy 27m = 3 x 10 51 erg, the duration AT = 30sec at the lOkpc, and the lower one 
to a GRB with 27m = 3 x 10 51 erg , AT = lOsec at the IMpc. The characteristic amplitude for the slope waveform with 
Ahm — 4\nvy/3 2 /r is used for the low-frequency region, / < 1/AT. The amplitude at / > 1/AT (plotted by the bold dashed 
line) depends on the structure of the GRB jet. The SNR can be obtained by equation 14.31 . 



The gravitational wave memory does not depend on the form of energy. For example, even if the jet starts out as 
a pure Poynting flux jet [TtI , the memory is observed in the similar way. One may also think that the memory is 
lost if the jet is decelerated, for example, at the internal shock, the external shock and the termination shock inside 
the massive star. However at the internal and external shock almost all kinetic energy of the jet goes into the photon 
emission whose memory also remains. (In equation (I2.1|) . 7m should be interpreted as the photon energy, and the 
photon speed (3 must be set to 1.) Thus the amplitude of the gravitational wave memory does not change so much. 
At the termination shock between the jet and the massive stellar envelope, the kinetic energy may be converted to 
the neutrino emission |49l l5l| . Since the neutrino also leave the memory, the gravitational wave memory may not be 
reduced so much as long as the cooling is not severe in the jet head 
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From equation 12.21 . we find that the ratio between the amplitude of the gravitational wave due to the forward jet and 

the one due to the counter jet is (1 + (5 cos9 v ) / '(1 — f3cos8 v ). Since we consider the case of 9 V < 0.25 in this paper, this 

ratio is much larger than unity. Therefore we can neglect the effect of the counter jet. 

Note that f™ x e l27rft H(t)dt = (l/2)[5(/) + Pv(i/7r/)] where Pv denotes the principal value 

The LIGO II noise curve is taken from We adopt S h (f) = 1.22 x 10" 51 /" 4 + 2.11 X 10" 41 + 1.22 x 10" 37 / 2 Hz -1 for 
LISA gjj, and S h {f ) = 4.5 x 10" 51 /" 4 + 4.5 x 10" 45 / 2 Hz" 1 for DECIGO/BBO 0, where / is in unit of Hz. 



